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Abstract We investigate polymer partitioning from polymer mixtures into nanometer size cav- 
ities by formulating an equation of state for a binary polymer mixture assuming that only one 
(smaller) of the two polymer components can penetrate the cavity. Deriving the partitioning equi- 
librium equations and solving them numerically allows us to introduce the concept of "polymers- 
pushing-polymers" for the action of non-penetrating polymers on the partitioning of the penetrating 
polymers. Polymer partitioning into a pore even within a very simple model of a binary polymer 
mixture is shown to depend in a complicated way on the composition of the polymer mixture and/or 
the pore-penetration penalty. This can lead to enhanced as well as diminished partitioning, due to 
two separate energy scales that we analyse in detail. 



I. INTRODUCTION 



First analyzed within the framework of size exclusion chromatography [l[ passive partitioning of polymers into 
nanoscale cavities has gained much broader relevance ■ Active forcing of polymers into nanosize pores underlies 
some applications of the "osmotic stress" technique, devised originally to probe inter- and intramolecular forces [|| but 
also used successfully to probe changes in the size of nano-cavities of some proteins [i[ . Size-dependent partitioning 
of water-soluble polymers is particularly important for the sizing and probing of water-permeable ion channels 0-01 
and for single- molecule mass spectroscopy [8-ll0l|. In addition, forcing polymers into nanosize tubes is relevant to 
recent studies of controlled ejection of viral genomes from capsids where osmotic stress is applied to push a 
DNA molecule into or to allow it out of the viral capsid by the action of PEG (poly-ethylene-glycol) dissolved in 
the ambient solution (T^]. In this last case, the identities of the polymer forced into the nano-cavity (DNA) and the 
polymer pushing it (PEG) differ. 





FIG. 1: A polymer solution composed of two polymer types (black and grey) in equilibrium with a pore. Only one component 
(black) of this binary polymer mixture is allowed to enter the pore, with a energy penalty A/. We assume throughout that 
the two polymer species differ only in length. 



We analyze a variation on the problem. Here the polymer solution is heterogeneous, composed of various sizes of 
the (same type of) polymers, for which the accessibility of the pore varies depending on their size. We thus assume 
that the polymer solution is polydispersc and that of the different size polymer chains only a single type can enter 
the pore (Fig. [T]). We show that in this case the polymers that are too big, or are otherwise prevented from entering 
the pore, act to push the other polymer into the pore. 

Usually penetration of polymers into molecular nano-pores incurs an energy penalty. It depends also on the 
concentration of polymer chains outside the pore, assumed to be in equilibrium with the pore. Several approaches 
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have been pursued in the theoretical elucidation of this problem [l|, [l?], [H| • The main difference between previous 
approaches and the one advocated here is that we start with a consistent free energy of the polymer-solvent mixture 
that allows us to formulate the thermodynamic equilibrium for the pore-external solution system in terms of the 
chemical potential of the polymer species as well as the chemical potential (osmotic pressure) of the solvent. This is 
important. We want to generalize the formalism in such a way that we will be able to deal with a mixture of polymers, 
whose components are selectively allowed to enter the pore. We introduce a simple ansatz for the free energy of a 
mixture of (the same type of) polymers varying only in their sizes, consistent with our previous phenomenological fit 
to the equation of state of bulk uncharged polymers fl6j . Calculating the partition coefficient allows us to ascertain 
that indeed some polymers can push others to enter the pore in a kind of osmotic tug-of-war. 

We build our analysis around the most simple case of a PEG polymer solution and first review the properties 
of a monodisperse solution. We then derive its equation of state, consistent with empirical fits over a wide range 
of concentrations (l6j . Then we generalize the equation of state of a monodisperse polymer solution to a simple 
polydisperse case, assumed to be composed of only two sizes of the same type of polymers: small and big PEGs, 
or sPEG and bPEG, respectively. The third explicit component of the solution is water. We will first analyze the 
osmotic pressure of this solution for various amounts of sPEG added to the background of bPEG. Then we assume 
that this solution is in equilibrium with a pore that can be penetrated by the sPEG at a finite free energy price but 
is impenetrable to bPEG. We calculate the corresponding partition coefficient and assess the pushing forces exerted 
by the external polymer solution via its osmotic pressure. 



II. FREE ENERGY OF A BINARY POLYMER MIXTURE 



We start with the free energy AF((f> p ) of a monodisperse PEG polymer solution as a function of the monomer fraction 
(ftp of polymers derived in [13j, whose range of validity is restricted to long polymer chains such that <\> p > <j> p *; <\> p * 
corresponds to the overlap concentration defining the onset of the semi dilute region, with monomer fraction low 
enough that higher order terms, above the second virial term, need not be taken into account. This form of the free 
energy and the consequent osmotic pressure is completely consistent with the phenomenological form of the PEG 
equation of state demonstrated in |16| . 

Under the stated constraintshe the free energy of this monodisperse polymer solution assumes the form [l3j 



AF(0 P ) 
k B T 
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where k B T is the thermal energy. The parameter a in the above free energy expression was evaluated in [1J 
a = 1.87; \ is the Flory-Huggins parameter and <p p is the monomer fraction of the polymer N p segments long 



n p N p _ n p N p 
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where n p is the number of polymer monomers, n w is the number of water molecules so that no is the total number of 
monomers plus water molecules. 

We now use this form of the free energy for a different problem, a binary polymer mixture. We assume that each 
of the two components separately conforms to the limits of validity of the above formula. If the mixture is composed 
of n s molecules of the small polymer N s monomers long, and n b molecules of the big polymer N b monomers long, in 
an aqueous solvent of n w water molecules, than the free energy of this polymer mixture is given by 

b - = n w hi (j> w + n s ln(j> s + n b \n(f) b + xn (4>s + <Pb) - ^n ((/) s + (f) b f + 
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a direct generalization of Eq. [TJ The definition of the monomer fractions is 

A nw A UsNs A UbNb l AT J AT ft\ 

9w = — <Ps = <Pb = with n = n w + n s N s + n b N b , (3) 

n n n 

so that (/> w + 4> s + <fib = 1- is the total number of molecules, i.e., the number of water molecules and all the monomers; 
the monomer fractions are expressed for long (b) or short (s) polymer molecules. 
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FIG. 2: (Color online) Mixture of PEG400 (N 3 ~ 9) and PEG 3500 (N b ~ 40). Plot of the log of osmotic pressure Eq. [5] as 
a function of the total polymer monomer fraction (log-linear plot). Four different values of <f> s are presented, viz. 4> s — (full 
line), O.f (O), 0.2 (□) and 0.3 (O). 



III. OSMOTIC PRESSURE 



The equation of state of the binary polymer mixture is its osmotic pressure as a function of the monomer fractions 
of the polymer types. We calculate this osmotic pressure H(<fi s ,<f)b) by first evaluating the chemical potential of water, 
for given n s and n&, 

^.™^»(**\ (4) 
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where V is the molecular volume of the solvent. The osmotic pressure of the solution in dimensionless form, H(4> s , <fib), 
is then obtained as 

fiP ^ vti{4> s m) 



k B T 

= -l n ^+^-l + ^ + A_i(i_^ + |a(i-x) 3/4 (l-^) 9/4 . (5) 

Remarkably, to the lowest order in <j) s and </>& this dimensionless osmotic pressure is a sum of the (ideal) van't HofF 
and the des Cloizeaux terms. In the limit of a single type of polymer, i.e. <f> s — or <pb = 0, this decomposition is 
consistent with the formula that was used to fit the equation of state experimental data in Refs. [IH, [l6j]. The only 

difference is in the notation: a — > a (| — x) 3 ^ 4 - Introducing a = a (| — x)"^ 4 one can deduce that a = 0.49 for 
PEG Hfj. 

Fig. [2] shows the osmotic pressure Tl(<f)s, <j>b) of a binary mixture of two PEG polymers of different sizes, PEG400 
with N s ~ 9 and PEG 3500 with Nb ~ 40, obtained from Eq. [5] For water V = 1.0 ml/g; we will assume the same 
value of molecular volume for PEG [2(j. We compare the osmotic pressure for 10, 20 and 30 % of PEG 400 on top of 
the background PEG 3500 polymer solution. As seen in Fig. [21 the osmotic pressure is higher for mixtures with the 
relative change decreasing with It must be noted that the effect of disparity in the molecular weights of the small 
and the big polymer, i.e., in N s and Nb, appears as the ratio of the total number of big molecules vs. the number of 
the small molecules, as is evident from Eq. [5] 



IV. CHEMICAL POTENTIALS OF POLYMERS 



For a mixed PEG solution in equilibrium with a pore that is permeable only to one component of the PEG mixture 
(e.g. sPEG) but impermeable to the other one (bPEG), the short PEG chain has to be in chemical equilibrium 
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between the solution and the pore [l[ . This signifies that the chemical potential of sPEG chain inside the pore and in 
the external solution must be the same. Thus it makes sense to define the partition coefficient of the PEG type that 
can penetrate the pore as a ratio between the concentration inside the pore and outside in the external solution. 

We first derive the chemical potential for the two PEG components of the solution as well as for water. For polymer 
s or b 

Vs.b = tt— [-, — ™ ) = In s .t, + 1 - (j> s b - (j> w N s h - (j>b, s -r^- - 
on S:b \k B T J N byS 

+(X - \)N s .b + \N s<b <j>l - f &N s , b {<t> s + h ) 9/4 + ^&N s . b {<p s + fof' 4 , (6) 

where, again, all the other variables are kept fixed in taking the derivatives. Using the expression for the chemical 
potential of water, Eq. |U we derive 

Hs,b-N s , b n. w = ln<£ s , 6 + l- N Stb {\n(j) w + l~(j) w ) + (x-l)N s .b + N s .bj&{(l) s +(j)b) 5/4 ■ (7) 

On the l.h.s. of the above equation we have an excess chemical potential of the polymer compared with the same 
number of solvent molecules as there are monomers. It is easy to show that the chemical potentials derived above 
satisfy the Gibbs-Duhem relation. 



V. PARTITION COEFFICIENT 



We will now calculate the partition coefficient for small polymers in solution assuming that they can penetrate the 
pore, schematically depicted in Fig. [TJ 

Assume first a polymer solution composed of only small (s) polymers and that sPEG can enter the pore with a free 
energy penalty A/. For the time being we surmise that this free energy difference does not depend on the state of the 
polymer in the pore and is independent of the polymer concentration inside the pore, depending only on the single 
chain confinement free energy. For the magnitude of A/, we use the standard argument [3, 0| based on the fact that 
the polymer needs to be squeezed into a pore whose diameter is smaller than the "natural" Flory radius of the chain. 
With this in mind wc obtain [ll, 0| 

Af = N s ^ = f (R)N s , (8) 

where a is the monomer size and where we assume that the length of the polymer that penetrates the pore is equal 
to the full length of polymer, N s . In what follows we will simply take the value of A/ to be independent of 4>s,b and 
do not specify its exact form. The latter indeed depends not only on the pore radius but also on the polymer - pore 
surface interaction energy, for details see Ref. . 

Then if / stands for the inside the pore and O for the outside, the chemical equilibrium is established when 

/i a (J) + A/ = /i s (0). (9) 

Also, since the solvent is in chemical equilibrium too, we need an additional equation 

M-0 = /MO). (io) 

In principle at least, for the solvent too we could add a separate pore penetration energy. For now we do not explore 
this scenario. The above two equations can then be rewritten in an equivalent form 

H„(T) - N sf i w (I) +Af = (x s (0) - N 8 fi w (0), (11) 

where we took into account the assumption that all the monomers and water molecules arc of the same size. Taking 
into account Eq. [6] 

ln0 s (J) - N s (ln^ w (I) + 1 - <M-0) + N s ^a (&,(/) + h (/)) 5/4 + A/ = 

in0 s (o) - A s (in <mo) + 1 - mo)) + N s ^a (UO) + MO)) 5/4 , (12) 
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FIG. 3: Plot of the partition coefficient Eq. [14] for a single species of polymer, PEG400 with N 3 = 9, as a function of the 
polymer monomer fraction <f) s {0) in the external solution. Cases with four different penetration free energies are presented, 
A/ = (full line), 2 (O), 4 (□), 9 (O) in dimensionless units. 

where we have removed all the irrelevant constants. Equation takes into account the chemical equilibrium of the 
pore-external solution system in terms of the polymer as well the solvent chemical potentials. Since, by assumption, 
the big polymer can not penetrate the pore we obviously have 4>b{I) = 0. Introducing now the partition coefficient 

Mil 

p Moy 

Eq. [TJ] assumes the final form 

lnp + A/ = N - ( ln (i -uo)-Zo)) + ipM0) ' UO) ~ M ° ))+ 

[(MO) + MO)) 5/4 - (P0-(O)) 6/4 )) . (is) 

The solution of this equation gives us p = p((f> s (0), A/) when solved numerically. Usually A/ is assumed to be a 
linear function of N s but this is irrelevant here. 

VI. PARTITION COEFFICIENT: SINGLE POLYMER TYPE 

Consider first that we have only small polymer, e.g. PEG400, in the pore as well as in the external solution, i.e. 
4>b(0) = 0. In this case, Eq. [13] reduces to 

InP + A/ = N s (in ^'^q)) + <* ~ ^ { ° ] + ¥ ^ ^) ^ [ ° f ") ' (W) 
Fig. [3] shows numerical solution assuming that a = 0.49 for PEG400. We evaluate the partition coefficient 

p = p(MO),Af) 

and plot it for different Af as a function of 4> s (0) in the external solution. 

The general features of the polymer partitioning in this case are well known: the partition coefficient is a mono- 
tonically increasing function of the monomer fraction of the polymer in the external solution, <p s {0), and strictly 
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p < 1, i.e. the monomer fraction inside the pore is always smaller than outside. The additional free energy penalty, 
A/, from polymer-pore interaction, makes it clearly more difficult for the polymer to enter the pore. Consequently 
higher values of the monomer fraction in the external solution are needed to reach the same partition coefficient. Not 
only that, as the penetration free energy penalty grows, the partition coefficient becomes a steeper function of the 
monomer fraction in the external solution, with a narrower interval of saturation behavior. 

We now evaluate explicitly two limiting cases for 4> s (0) <C 1 and </> s (0) — > 1. In the first case <fi s (I) = p<fi s (0) < 
4> s (0) so that obviously p <C 1, we obtain to the lowest order 

hip + A/ = N s (|« (1 -p 5 / 4 ) MO) 5/4 + OlMO)} 2 ) s N s \a <f> s (0) 5/4 - (15) 

and therefore 

\np(<j> s (0)) ~ - A/ + N s j5t MO) 5/4 - (16) 

The partition coefficient is thus increased from its ideal value given by lnp(0 s (O)) = lnpo = — A/. In the opposite 
limit of tp s (0),p — > 1, within the range of \p — 1| -C A/, 

]np{UO)) * -A/(JV.fS UOf /4 Y\ (17) 

i.e., the partition coefficient is decreased from its ideal value. These two limiting cases can be well discerned also 
from the complete numerical solution of Eq. [14] (Fig. |3j). It is clear from this figure that the range of validity of the 
limiting expression Eq. 1171 depends on the value of A/. 

VII. PARTITION COEFFICIENT: TWO POLYMER TYPES 

We now return to the original problem of a polymer mixture, in the simplest case composed of two types of polymer, 
b and s, of which only the s can penetrate the pore. In this case, we need to solve the full equation Eq. 1131 The 
solution 

p = p(<t> a (0),MO),&f), (18) 

is a function of the monomer fraction of the small and the big polymer in the external solution as well as of the 
free energy penalty for entering the pore. We plot the solution p = p((p s (0),<pb(0)) for different values of A/ as a 
parameter, which corresponds to the variation of the penetration free energy with the pore size. As an elucidating 
example, we again have for "s" polymer PEG400 with N s = 9 and for "b" polymer, PEG 3500 with N b = 40. We 
study separately the dependence of the partition coefficient on <t>b{0) with 4> s (0) fixed, and vice versa. Both cases 
differ fundamentally from the single polymer type pore partitioning. 

Fig. 2] shows the dependence of the partition coefficient on the composition of the external solution, <p s (0), <f>b(0), 
as well as the free energy penalty A/. Consider first the dependence of the partition coefficient on (f>b(0) at fixed 
4> s (0), Fig. QJa), (b). We see that at every set 4> s (0) the addition of the non-penetrating polymer "b" to the external 
solution monotonically increases the partition coefficient of the penetrating one, (Fig. Eta)). We can therefore 
conclude that the "b" polymer always pushes the "s" polymer into the pore. This is the physical meaning of the 
"polymers pushing polymers" concept. In addition, for every composition, i.e., <p s (0),<pb(0), there exists a "critical 
value" of the penetration free energy penalty at which the monotonically decreasing partition coefficient crosses the 
value p = 1, see Fig. [5{a), defined as Afi(<p s (0), <f)b{0)) = Af(p = 1), where we indicated explicitly that it depends 
on the composition of the external polymer solution. Depending on the composition of the external solution and for 
A/ < A/i, we can have enhanced partitioning of the penetrating polymer with p > 1. This is very different from 
the case of a single polymer type where always p < 1. Note also that the derivative of the partition coefficient with 
respect to 4>b{0) is always positive, for any composition and any penetration penalty. The partition coefficient can 
be analyzed alternatively also through its dependence on the total osmotic pressure of the external solution, obtained 
from Eq. [5l (Fig. Hfb)). Clearly for fixed 4> s (0) the polymer "s" is always pushed into the pore, Fig. 0Jb), as its 
partition coefficient is an increasing function of the total external solution osmotic pressure. 

The converse case of the dependence of the partition coefficient on <fi s (0) with constant (j>b(0) is more complicated. 
First of all, in this case too we discern, (Fig. SJc)), that polymer "s" can have enhanced partitioning, i.e., its partition 
coefficient can be larger than one. The corresponding "critical value" of the penetration free energy penalty is the 
same as in the previous case, Afi(4> s (0),4>b(0)) = Af(p = 1) and is an increasing function of both 4> s {0) and 
4>b{0), but is not symmetric in the two arguments, as can be clearly seen from Fig. IHta). There is, however, an 
additional feature in the 4> s (0) variation of the partition coefficient when 4>b(0) is kept fixed. From Fig. HJc), we 
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FIG. 4: (Color online) (a) Partition coefficient Eq. \W\ for different values of 4>b(0) and different values of the additional 
polymer-pore interaction free energy A/, at a fixed value of <f> a (0) in the external solution. <f> s (0) = 0.05 and A/ = 0, 4, 9 (O, 
□ , O). (j> s (0) = 0.1 and A/ = 0,4,9 (•, ■. ♦). (j} s (0) = 0.2 and A/ = 0,4,9 (®, EB, M). The delimiting value p = 1 is shown 
as a black line. 

(b) Partition coefficient for different values of 4>b{0) and the same fixed value of 4> s {0) in external solution, as a function of 
total dimensionless osmotic pressure tl((j> s (0), <f> b (0)). Numerical results for three different values of A/ and three different 
values of 4>s{Q) are presented. 

(c) Partition coefficient p((f> a (0), cj> b (0), A/), Eq. 1161 for different values of PEG400 (N b = 9) monomer fractions 4>a{0) and 
different values of A/, at a fixed value of PEG3500 (N b = 40) <j> b in the external solution. <j>b(0) = 0.2 and A/ = 0, 4, 9 (O, □, 
O). <t> b (0) = 0.1 and A/ = 0,4,9 (•, ■, ♦). tf> b (0) = 0.05 and A/ = 0,4,9 (©, ffl, H). The delimiting value p = 1 is again 
shown as a black line. 

(d) Partition coefficient for different values of 4> a (0) at the same fixed values of 4> b (0) in external solution, as a function of 
the total dimensionless osmotic pressure fl((f> s (0), </> b (0)). Numerical results for three different values of A/ are presented and 
three different values of <j> b (0). 
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FIG. 5: (a) Partition coefficient Eq. [Til for PEG400 (N s = 9) as a function of A/, for cj> a (0) = 0.2, <f> b (0) = 0.5 (O), 
<?!> s (0) = 0.2, tj> b (0) = 0.3 (□), 4> s {0) = 0.2, b (O) = 0.2 (A) and <j> s (0) = 0.2, tj> b (0) = 0.05 (O). Partition coefficient thus 
appears as a monotonically decreasing function of the pore penetration energy A/, (b) Derivative of the partition coefficient 
a<MO) \ <t>b(0) f rom Eq. [19] as a function of A/, for different external solution compositions: <fis(0) = 0.2, <j> b (0) = (A), 
<t> 3 {0) = 0.3, <f>b(0) = 0.1 (□), <f>s(0) = 0.4, (j)b(O) = 0.2 (O). The derivative changes sign at a "critical point" of the pore 
penetration free energy, A/ c ~ 1 in the second case and A/ c ~ 2 in the third case. For large enough A/ the derivative is 
positive and limits to zero. 
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FIG. 6: (a) "Critical value" of the penetration free energy penalty defined as A/i(0 s (O), </> b (0)) = A/(p = 1) for each particular 
composition of the binary mixture (j> s (0),4> b (Q). As long as A/ < A/i there is enhanced partitioning or superpartitioning of 
the penetrating polymer with p > 1. Afi(cj> s (0), 4> b (0)) is a monotonic function of the composition, (b) "Critical value" of the 
penetration free energy penalty defined as Afo((j>s(0),<f> b (0)) = A/(p' = 0). For A/ < A/o we have p = g ^ ) ^ ^ an< ^ 
opposite for A/ > A/o, thus leading to augmented and depleted partitioning regions in the parameter space. Afo(4> s (0), (j> b (0)) 
is obviously a non-monotonic function of the composition. 



discern that for any non-zero value of c/>b(0), there is a range of Af values where the rate of variation of the partition 
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coefficient with respect to <f) s (0), i.e., p' = g^jp) > cnan S es sign. Apart from possible enhanced partitioning we 
can thus also have an augmented, p' > 0, and a depleted rate of partitioning with p' < 0. In complete analogy 
with the behavior of the partition coefficient itself, we can define a "critical value" of the penetration free energy as 
Afo(4> s (0), <j>b(0)) = A/(p' = 0), delimiting the augmented and the depleted partitioning regions, sec Fig. [5jb) and 
Fig. Mb). 

Fig. EJb) also shows that Afo((j> s (0), 4>b{0)) is not a monotonic function of its arguments. If we fix 4>b(0) and follow 
the variation of A/o with <f> s (0), we note that it first increases then decreases, while it appears to be strictly monotonic 
in the converse case of fixed (f> s (0) and variable 4>b(0). The non-monotonic variation of A/o(0 s (O), 4>b(0)) indicates 
a complicated interplay of various parts of the free energy on polymer penetration of the pore, specifically of the 
polymer mixing entropy and configurational fluctuations in the external binary polymer solution and the penetration 
penalty for entering the pore by only one type of the polymers. There is no analogous effect in the case of a single 
polymer type inside and outside. 

Since the depleted partitioning of polymer "s" form the pore, as defined above, is a counterintuitive phenomenon, 
it needs some further elucidation. For every composition of the system, that is for every fixed pair (f) s (0),4>b{0), 
there is thus a corresponding Af ((f> s (O), <pb(0)), where the rate of variation of the partition coefficient with respect 
to 4» s (0), changes sign. In fact one can derive the following simple relationship 

d\ogP , _ MO) dp _ N s (xs(Q)-Xs(I)) ncn 

d\ogMOy MO) P dUOy MO) l + N sXs (I) ' [ ! 

where we introduced the osmotic compressibilities of the penetrating polymer as 

^ = ^mry (20) 

Here 4> s {* = O, I) stands for the monomer fraction outside or inside of the pore, while II = H((j> s (0, I), <frb{0, 1)) is the 
corresponding dimensionless osmotic pressure, with both osmotic compressibilities depending implicitly on A/. From 
Eq. [T2]it follows that the rate of variation of the partition coefficient with respect to 4> s (0) depends on the difference of 
the osmotic compressibilities of the penetrating polymer outside and inside the pore. The transition from augmented 
to depleted partitioning is thus connected with the osmotic compressibility of the exchangeable polymer in the two 
environments and one could say that A/o(0 s (O), <f>b(0)) is defined also by the equality of the two comprcssibilituies, 
Xs{0) and Xs{I)- 

The behavior of the polymer partitioning in the case of variable <p s (O) is therefore governed by two energy scales 
A/o and A/i determining the rate of variation and the magnitude of the partition coefficient p(4> s (0), <j>b{0), A/). 
It is quite unexpected that the mere introduction of an additional non-penetrating polymer to the external solution 
would have such fundamental and non-trivial repercussions for the partitioning of the other polymer. 



VIII. DISCUSSION AND CONCLUSIONS 



We present a thermodynamic equilibrium analysis of pore penetration by polymer chains in a binary polymer 
mixture, where only one of the components is allowed to penetrate the pore with a set free energy penalty. 

First we introduced an equation of state, i.e., the dependence of osmotic pressure on the monomer fraction of the 
polymers, Eq. [5J consistent with empirical fits obtained for the most complete set of experimental data available for 
PEG in aqueous solvent as well as PAMS in toluene The proposed equation of state captures all the features of 
the empirical fits. We generalized the equation of state for a single type of polymer to the case of a simple binary 
polymer mixture, composed of two types of polymer chains, "s" and "b", of which only the former can enter the pore 
with a finite free energy price. We then used this equation of state in order to formulate the equilibrium distribution 
of the short polymer between the external solution and the pore, obtaining the partition coefficient of the short chains 
in terms of the external polymer solution composition. Our analysis reduces to the known cases when there is only a 
single type of polymer in solution. 

The partitioning of a single component of the polymer solution into a pore can be described with the concept of 
"polymers pushing polymers": the non-penetrating polymer chains "b" contribute an osmotic push for an enhanced 
partitioning of the " s" polymer chains into the pore. We find that the larger the concentration of the non-penetrating 
" b" polymer in the external solution, the larger the concentration of penetrating polymer " s" that will be pushed into 
the pore at any value of the pore penetration energy. However, this does not imply that the partitioning coefficient, 
i.e., the ratio of the "s" concentration inside and outside the pore, is always a monotonic function of the external 
concentrations of both polymer types. 
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Analysing the general properties of the partition coefficient when the "s" and the "b" components of the external 
polymer mixture are allowed to vary independently, we find that the general features of the pore penetration are 
governed by two different penetration free energy scales, A/i and A/o. The former differentiates between enhanced 
partitioning, p > 1, and ordinary partitioning p < 1 dependent on the osmotic pushing strength of the external 
polymer mixture and indeed embodies the concept of " polymers pushing polymers" . The other penetration free 
energy scale, A/o, that differentiates between the rates of variation of the partitioning coefficient of the penetrating 
"s" polymer leading to an augmented, p' > 0, and a depleted rate of partitioning, p' < 0. While A/i(</> s (0), (f>b(0)) 
is a monotonic function of its arguments, A/o(0 s (O), 4>b(0)) is not. The latter type of behaviour is understood as 
stemming from the interplay between the polymer mixing entropy and/or configurational fluctuations in the external 
binary polymer solution, and the pore penetration penalty. 

Both of these free energy scales depend on the composition of the external solution and for a fixed A/ one should 
be able to control the partitioning of a polymer into a pore, or in general into any nano-cavity, by varying the 
composition of the external polymer mixture. This principle should have applications in many areas of nano-science 
dealing with partitioning of polymer chains into small enclosures. Exp eriments are under way to test at least some of 
the conclusions reached on the basis of the theory presented here [13]. 




While we have analyzed only the most simple case of a binary polymer mixture, one can in principle generalize this 
theory to more complex polydispersity. This would entail the introduction of a size distribution in the equation of 
state as well as an energy size distribution for the pore penetration penalty. One could formulate the corresponding 
expressions for equilibrium partitioning in a parallel fashion to what we derived for the binary mixture. This would 
open the possibility to address new problems in the theory of confined polymer solutions. Furthermore it is easy to 
envision cases where the polymer mixture would be actually composed of polymers of very different nature such as, 
e.g., PEG and proteins, or PEG and sugars. It is reasonable to expect that a similar analysis would apply to those 
cases too. 
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